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I. INTRODUCTION 



Recently, the Riccati operator equation (RE) has attached an attention to mathematical 
physicists (see e.g, Refs.^''^). As it has been argued, the solution of the equation can be used 
to obtain reduced time evolution pt of a two-level open quantum system'''. This, so-called 
the reduced dynamics, plays a central role both in the quantum information theory and the 
theory of quantum information processing. It is well-established that pt = ^t{po), where 

is the trace preserving (TP) and completely positive (CP) map ' . Such TP-CP maps, 
known as channels, are commonly used to simulate real quantum systems and to design 
quantum algorithms''*'. 

To understand a connection between a given TP-CP map $4 and a solution of the RE, let 
us begin by reviewing some basics facts concerning open quantum systems. For this purpose 
we consider the system-environment Hamiltonian in the following form 

H = ffQ®Is + lQ®i/E + Hint, (1) 

where Hq and He represent the Hamiltonian of the qubit and the environment, respec- 
tively. Hint specifies the interaction between the systems, whereas Ig and I^; are the identity 
operators on corresponding Hilbert spaces and He of the qubit and its environment, 
respectively. H acts on the Hilbert space "Htot = 'He and it admits the block operator 
matrix representation'': 

V 

H = 

H 

with respect to the decomposition "Htot = "He ©"He- The diagonal entries H± : V{H^) — )■ He 
are densely defined and self-adjoint unbounded operators. We have assumed that "DlV) = 
V{V'') = He, which means that V, and thus V'' as well, is bounded. 
The time reduced evolution of an open system is given by 

$i(po) = TrE[Uf^(po)Ul], (3) 

where Ut = exp(— zHt) is the time evolution operator of the total system. The map 
assigns to each initial state po a single state \l/(po) of the total system. This assignment map 
must be chosen properly so that $i can be well-defined^". In this paper we assume that no 
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on ViH) = V{H^) © V{H_ 



(2) 



correlations between the systems are initially present^^ and thus we take \l/(po) = Po ® 
for some initial state of the environment u. 

The linear map Tt^ denotes the so-called partial trace: 



Tr, 



Mn 


M12 




M21 


M22 





G M2 (C) , where Mij G T ("He) ■ 



(4) 



TrMii TrMi2 
JrM2i TrM22 

In the above description, Tr refers to the usual trace on "He, T('He) denotes the Banach 
space of trace class operators with the trace norm: \\A\\i = Tr(vAAJ), whereas M2(C) is 
the Banach space of 2 x 2 complex matrices. Form (4) one can see that the partial trace in 
Eq. (3) transforms states of the total system (block operator matrices, square bracket) into 
states of the open system (complex matrices, round bracket). It is worth noting that the 
partial trace TrE, unlike the trace Tr, is not cyclic, namely TrE(AB) 7^ TrE(BA). 

In general, the formula (3) is far less useful, than its simplicity might indicate. Indeed, 
to trace out the state Uj\l/(po)Uj over the environment degrees of freedom one needs to find 
a decomposition of Ut with respect to the decomposition of Hilbert space "Htot = ^e © ^e- 
This problem is fundamental in quantum mechanics, as it requires the diagonalization of H. 

It is well-established (see e.g^^ and Refs. therein) that H can be block-diagonalized 
provided there exist a bounded solution X (with Ran{X)\x>(H+) C Vi^Hj)) of the Riccati 
equation: 



XVX + XH+- H_X -V^ = Q on V{H+). 
To be more specific, if X solves (5) the following equality holds true: 



S^HS 



Z+ 

z_ 



where S 



(5) 



(6) 



X Ie 

and Z+ = H+ + VX, Z_ = H_ - V^X^ with V{Z^) = V{H^). Moreover 8"^ exists as a 
bounded operator'. 

Using decomposition (6) one can write in an explicit matrix form: 
















Ut = exp 1 

















z_ 







-iZ+t 











iZ-t 



(7) 



Among all difficulties concerning the RE, we will focus on two problems described bellow. 
The first one^ arises when the solution of the RE is antilinear. In this case, the formula (7) 



cannot be applied, since S is neither an antilinear nor a linear operator. Therefore, we are 
interested in asking how can we determine the time evolution of the total system. 

The RE is an operator equation and even for a simple system it may be problematic to 
find the solution. However, there are cases in which the RE can be simplified considerably^ . 
For instance, in the case when V = alE = a, where a G M; then the RE reads 

aX'^ + XH+- H_X -a = on V{H+). (8) 

In the instance when the operators iJ± are self-adjoint and commute [e.g, [H^, H_] = 0), it 
is reasonable to assume, that the solution has the following form X = or X = g{H_). 

Then, the functions / and g read / = h\^(^H+) and g = /i|cr{H_), where h satisfies the quadratic 
equation: a/i(A)^ + 2h{X)X — a = 0. At this point the second issue can be addressed: does 
every solution of (8) is an analytical function of or H_7 

The primary goal of the presented work is to solve this two worthwhile issues. First, we 
show that there is a simple way of computing the evolution operator with the use of 
a formula similar to the one provided by the Eq. (7). Next, we exemplify that not every 
solution of (8) is an analytical function of or H_. 

II. ANTILINEAR SOLUTION 

In this section we derive a simple method of computing in the case when the solution 
of the RE is antilinear. To end this let us assume that r is a bounded, antilinear (z.e, 
T^zlip)) = z*t\'iIj)) operator which solves Eq. (5) and let Sr = where S denotes matrix 
from Eq. (6) with X = t. Although, Sr is neither a linear nor an antilinear in this case, it 
has an intersecting property: 

Sr{z\^)) =^{z)Sr\ip) +iQ{z)S_r\^) (uotc S_,S, = 1^ = S,S_,) , (9) 

which reduces to Srizli/j)) = zSrli'), for z G M. Therefore, S,- can be treated as a linear 
operator when it acts on some subspace over a field of the real numbers M. In other 
words, the restricted operator St-i-^^ is linear. 

As a result, if H = ST-HdS_T- then for an analytical function / : o"(H) — )• M we have 
/(ST-HdS_T-) = ST-/(Hd)S_^, thus in particular we obtain 

exp(— iHt) = St- cos(Hdt)S_T- — iS,- sin(Hdt)S_T-. (10) 
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Although the right side of the Eq. (10) is not equal to S,- exp(— 'iHdt)S„T-, yet the evolu- 
tion operator can still easily be managed by separately transforming cos(Ht) and sin(Ht) 
functions, and then combing the results by using the Euler formula. 



III. NON-ANALYTICAL SOLUTION OF THE RICCATI EQUATION 

In this section we give an example that the solution of (8) is not necessarily an analytical 
function of or H_. To achieve this, we consider the RE (8) with = ±H, where 
H : Vi^H) — )• "He is assumed to be densely defined and self-adjoin operator. With the use of 
these assumptions the RE oiaV{H) can be written as 



If one assumes that the solution is in the form X = f{H), then the RE (11) simplifies on 
T){H) to the quadratic operator equation: 



It is not difficult to see that the function / is a solution of the quadratic equation: 



Since there is only one self-adjoint operator in Eq. (11), the assumption X = f{H) seems 
to be reasonable. The question concerning analytical solutions of RE in this case reduces to 
the question: are the Eqs. (11) and (12) equivalent? We will exemplify that the answer in 
general is negative, and we will prove it choosing H such that the solution of (11) does not 
solve (12). 

We construct H as follows. Let and a be the bosonic creation and annihilation op- 
erators, respectively. They are defined on a common domain Vi and obey the canonical 
commutation relation: [a, a"!"] = I, on some dense subset V2 of Vi. In terms of this operators 
we define the number operator = a^a, which is self-adjoint on V2. An explicit form of 
the sets is given by^'^ 



aX^ + XH + HX-a = 0. 



(11) 



aX^ + 2HX -a = 0. 



(12) 



a/(A)2 + 2A/(A)-a = 0, where \ea{H). 



(13) 




k 



1,2; 



(14) 
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where the vectors G 1^2 satisfy A^|0n) = '^l^n)- Moreover, on Vi the creation and 
annihilation operators can be defined exphcitly as 

oo oo 

a\ct>) = ^ v^(0„|0)|0„_i), at|0) = ^ v/^(0„|</>)|0„+i), 10) e I^i- (15) 

n=l n=0 

Finally, we take H = g*a + ga'' with T>{H) = Vi, where g E C 

It is not difficult to see that the solution of the RE (11) is given by the well-known bosonic 
parity operator 

oo 

P\ij) = J2^'''''{<f>n\mn), \^)eHE. (16) 

n=0 

Indeed, from (16) follows that P is both Hermitian and unitary, in particular = Ie- 
Therefore, in order to prove that P solves (11), it is sufficient to show that PHP = —H. 
The latter is obvious since PaP = —a and Pa) P = —a) . 

Note, formally P can be written as P = exp(i7rA^), however, unlike it is everywhere 
defined, thus it is bounded as well (||P|| = 1). Of course, P is not an analytical function of 
H and it does not solve Eq. (12) because PH ^ 0. 

IV. SUMMARY 

We have proposed a simple resolution of the recent problem concerning antilinear solution 
of the Riccati equation. Our result allows to take advantage of the result (7), in the case 
of an antilinear decomposition of the total Hamiltonian. The solution given in this paper is 
not as complicate as it is crafty. It relays on the very same exponential formula that is uses 
for a linear decomposition of the total Hamiltonian. However, the way it is used is different 
as we have explained in Sec. II. 

We have also argued that a simpliffcation of the RE, even under reasonable assumptions, 
can lead to the not equivalent equation. 
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